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Abstract. Observational determinations of the correlation be- 
tween the characteristic density <5 C and the virial mass M v ; r of 
dark halos constitute a critical test for models of hierarchical 
structure formation. Using the dynamical properties of dark 
halos reconstructed from the galaxy distributions in massive 
systems (groups/clusters) and the rotation curves in less mas- 
sive systems (dwarf, low surface brightness and spiral galax- 
ies) drawn from the literature, we confirm the existence of 
the 5 C -M v j r relation over a broad mass range from 10 10 M Q 
to 10 15 M Q , which is in gross consistency with the prediction 
of a flat cosmological model with £1m = 0.3 and = 0.7. 
It is pointed out that previous analyses based on the measure- 
ments of X-ray emitting gas and the hydrostatic equilibrium 
hypothesis, which claimed a shallower scale-free spectrum of 
n ~ — 1 for initial density fluctuations, may suffer from non- 
gravitational heating influence especially in low-mass systems. 
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1. Introduction 

In hierarchical clustering models, low-mass halos (e.g. galax- 
ies) are denser and have higher collapse redshift than massive 
halos (e.g. clusters), while the latter form by the gravitational 
aggregation of individual low-mass objects. For a halo density 
profile of the form Pdm(^) = ^cPcP( r / r s) where p c is the 
critical density of the universe, <5 C and r s are the characteris- 
tic density and length, respectively, the above scenario implies 
that the characteristic density 5 C of virialized dark halos should 
correlate with virial mass M v i r . Indeed, the existence of such 
a <5 C -M v i r relation has been confirmed by a number of high- 
resolution A^-body simulations in various cosmological models 
(SCDM, LCDM and OCDM) with different power spectra (e.g. 
Navarro et al. 1997; Salvador-Sole et al. 1998). 

On the observational side, several attempts have been re- 
cently made to test the <5 c -M V i r relation and hence, to constrain 
the scale-free power spectrum (P(k) oc A;™) for primordial den- 



Send offprint requests to: L.-M. Cheng 
Correspondence to: clm@class2.bao.ac.cn 



sity fluctuations on relevant scales. By fitting the X-ray sur- 
face brightness profiles of clusters predicted by the universal 
density profile (Navarro et al. 1995; NFW) to the ones for 63 
rich clusters observed with ROSAT, Wu & Xue (2000) found 
that 5 C oc M~^' 2 , indicating n i=s —0.7. This result was sub- 
sequently confirmed by Sato et al. (2000) based on an anal- 
ysis of the mass profiles of 83 ASCA X-ray objects over a 
wider mass range from M vu — 1O 12 M to 10 15 Af Q , which 
gives n = —1.2 ± 0.3 (90% confidence level). It appears that 
while the presence of good correlation between S c and M v [ r 
has been observationally established, the resulting power spec- 
tra are somewhat shallower than the value (n w —2) expected 
for typical CDM, scale-free cosmologies on cluster scales. Pos- 
sible reasons for this conflict have been outlined in Mahdavi 
(1999) and Wu & Xue (2000). Basically, the above observa- 
tional tests of the <5 C -M v i r relation based on X-ray observations 
of the intragroup/intracluster gas may suffer from the influ- 
ences of preheating, cooling flows, non-thermal pressure and 
hydrostatic equilibrium hypothesis. Among these, preheating 
by supernovae and AGNs in the early phase of structure for- 
mation may lead the intragroup/intracluster gas to extend out 
to larger radii (e.g. Kaiser 1991; David et al. 1991; Ponman 
et al. 1999; Bower et al. 2001), which can significantly affect 
our mass estimates from the X-ray observed surface brightness 
profiles especially for galaxies and groups when incorporated 
with the hydrostatic equilibrium hypothesis. 

Unlike intracluster gas, galaxies in groups and clusters are 
essentially unaffected by the presence of preheating, cooling 
flows and non-thermal pressure and can, therefore, be regarded 
as better tracers of underlying gravitational potentials. In less 
massive systems such as the dark matter dominated dwarf 
galaxies or low surface brightness galaxies, one is able to re- 
cover the dark matter profiles from the well-measured rotation 
curves (e.g. Carignan & Beaulieu 1989). Dynamical proper- 
ties of dark halos derived from the distributions of galaxies in 
groups/clusters and the rotation curves of dwarf, low surface 
brightness and spiral galaxies, which spans several decades in 
mass, have already been available in the literature. It thus be- 
comes possible and also timely to have a close examination of 
the <5 C -Af v i r relation motivated by hierarchical clustering mod- 
els. We attempt to fulfill the task in this paper, and compare our 
new determination with the results given by X-ray observations 
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and theoretical expectations. Throughout this paper, we assume 
H a = 100 km s _1 Mpc -1 and a flat cosmological model of 
Qm — 0.3 and Oa = 0.7. The shape of the <5 C -M vir relation is 
roughly unaffected by this choice (Wu & Xue 2000; Sato et al. 
2000). 

2. Sample 

We model the virialized dark halos by the NFW universal den- 
sity profile: pdmM = S c p c /[(r/r s )(l + r/r s ) 2 ], where the 
critical density of the universe reads p c = (3Hq/8itG)E 2 (z) 
and E 2 {z) = (1 + z) 2 {l + zfl M + [(1 + z)- 2 - l]ft A }. The 
virial mass M v - U is defined by M v - lr = 47rr^ ir A c p c /3, so that 
6 C = (A c /3){c 3 /[ln(l + c) - c/(l + c)]}, where A c is the 
overdensity of dark matter with respect to the critical value p c 
and can be approximated by A c = 178f2M(z) ' 45 (Eke et al. 
1998), and c = r v ; r /r s is called the concentration parameter. 

Assuming that galaxies in groups and clusters trace mass, 
we are able to fix the free parameter r s by straightforwardly fit- 
ting the observed surface number density profiles of galaxies to 
the projection of the NFW profile. In order to obtain 5 C , we also 
need to have an independent estimate of the virial radius r v ; r 
or the concentration parameter c. An approximate approach to 
the problem is to employ virial theorem: M v ; r = Sa 2 r v -„/G, 
which yields r v ; r = (6/ A c ) 1 ^ 2 [a p /HoE(z)], where cr p is the 
global line-of-sight velocity dispersion. However, this expres- 
sion holds true only for clusters with an isotropic distribution of 
galaxies, and their velocity dispersion can be reliably measured 
out to virial radii. Recall that the velocity dispersion of galaxies 
varies with radius if galaxies are assumed to trace mass (Cole 
& Lacey 1996). We have thus checked the accuracy of the es- 
timated virial radii using a sample of 158 nearby clusters of 
Girardi et al. (1998), in which the virial radii have been given 
by dynamical analysis, and found that there is good agreement 
between the two methods. This technique has been applied to 
an ensemble of rich systems of galaxies by CNOC (Carlberg 
et al. 1997), and poor systems of galaxies (PSG) by Mahdavi 
et al. (1999). Their averaged, best-fit parameters are summa- 
rized in Table 1 . For PSG we adopt the results of the sample D 
[PSG(D)] and eight groups with a p > 350 km s" 1 [PSG(H)] 
but exclude the result for the nine groups with <r p < 350 km 
s _1 due to the failure of the NFW fit. The third data set is taken 
from the ENACS cluster sample fitted by Adami et al. (1998). 
These authors used a similar technique but a pseudo NFW pro- 
file. We adopt their average value of r s = 0.26 Mpc among 
the 41 clusters whose velocity dispersions are observationally 
determined. The recent work based on a sample of 77 compos- 
ite ENACS clusters gives the similar result, r s = 0.318 Mpc 
(Adami et al. 2001). 

For less massive systems, the two free parameters in the 
NFW profile can be determined by utilizing the high-quality 
rotation curves, van den Bosch & Swaters (2001) have recently 
obtained the concentration parameters c 2 oo and the circular ve- 
locities V200 at r2oo for a sample of 20 dwarf galaxies by prop- 
erly subtracting the contributions of the thin gas disk and the 
thick stellar disk to the rotation curves. We take the mean val- 



ues of their best-fit C200 and V200 f° r the 15 dwarf galaxies 
(DWG) that satisfy c 2 oo > 1 and convert c 2 oo and r 2 oo into c 
and r vir in terms of the NFW profile and our definitions of c and 
r v ; r . The major uncertainty in their fittings arises probably from 
the constant mass-to-light ratio M/ L assumed for the stellar 
disk components. In order to demonstrate this uncertainty, we 
adopt both the result for M/L = l(M/L) Q (DWG1) and the 
one for M/L = (DWG0). Another sample we have used is 
the 9 low-luminosity disk galaxies (DIG) studied by Borriello 
& Salucci (2001). These authors have worked out the c and r s 
parameters based on the fitting of the rotation curves to the stel- 
lar disk plus dark halo models, in which the mass-to-light ratio 
M/L for the stellar disk is treated as a free parameter. However, 
they found that good fits in DIG sample are obtained only for 
unreasonably large virial velocities and masses. To overcome 
the inadequacy, they imposed an upper limit of 2 x 10 12 M© 
on M v ; r in their fitting. We use their average values of c and r s 
after the cosmological model correction is properly made. We 
have also performed the fitting of the rotation curves predicted 
by the NFW halo plus thin stellar disk model to a sample of 
30 spiral galaxies in the Ursa Major (UMA) cluster analyzed 
by Sanders & Verheijen (1999). However, the models are not 
well constrained for majority of galaxies due to the sparse data 
points of the rotation curve if the mass-to-light ratio of stellar 
components is treated as a free parameter and no upper limit is 
imposed on the virial masses of dark halos. We have only used 
the results of 13 galaxies that satisfy c > 1, and their best-fit 
parameters are listed in Table 1 . 

As a comparison, we turn to the best-fit values of r s and 
c from the X-ray measurements of 63 rich clusters observed 
with ROSAT (Wu & Xue 2000) and 83 objects observed with 
ASCA (Sato et al. 2000). While these two groups both assumed 
hydrostatic equilibrium for X-ray emitting gas, they used very 
different methods in the determinations of r s and c: Wu & Xue 
(2000) directly fitted the theoretically predicted X-ray surface 
brightness profiles Sx (r) from the NFW profile via isothermal- 
ity (Makino et al. 1998) to the ROSAT observed S x {r); Sato et 
al. (2000) first derived the mass distributions of groups/clusters 
from the X-ray observed surface brightness profiles character- 
ized by conventional (3 model and the temperature profiles. 
They then calculated <5 C and r s by fitting the NFW model to 
the mass profiles. For our analysis below, we use 3 and 4 bins 
according to virial mass for the ROSAT and ASCA samples, 
respectively (see Table 1). 

3. Analysis and results 

Fig. 1 shows the S c -M v i r relation over a broad mass range from 
10 10 M Q to 1O 15 M , in which the data points from X-ray emit- 
ting gas and those from distributions and rotation curves of 
galaxies are clearly marked. All the quoted error bars are 95% 
confidence limits. Also plotted in Fig.l is the corresponding 
c-M v i r relation for the purpose of easier comparison with nu- 
merical results (e.g. Navarro et al. 1997). It can be seen that the 
two data sets (galaxies and X-ray gas) are essentially consistent 
with each other on large mass scales M v ; r > 10 14 M Q , but dif- 
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Table 1. Sample 



sample 


No 


redshift 


<T p (kms x ) feT(keV) 


r a (Mpc) 


c 


M vir (1O 14 M ) 


5o (10 4 ) 


ref 


CNOC 

ENACS 

PSG(H) 


14 
41 

8 


0.313 ± 0.120 
0.073 ± 0.027 
0.025 ± 0.000 


972 ± 48 
617 ±38 
460 ± 38 


0.26 ± 0.03 
0.17 ± 0.05 


3.701!;" 

5-42l ;- 

6.00 ± 1.70 


11 76+ 1 ' 82 

' '"-1.65 
, 7 „+0.73 
•'■ ' U -0.65 

1 60+ 043 


0.29t ;f 8 

o,58+° ; 34 

7Q+ - 56 


1 

2 
3 


PSG(D) 


5 


0.024 ± 0.000 


392 ± 48 


0.117 ± 0.043 


6.70 ± 2.40 


0.99±°;£ 


91+ 98 

u - J -0.58 


3 


UMA 
DIG 


13 
9 


0.004 
~ 




0.022 ± 0.002 
0.023 ± 0.002 


9.71 ± 0.45 
7.37 ± 0.62 


8A6t° ° a 9 g x lO" 3 
5.61^ f 7 x 10- 3 


, q-,+0.48 
171 —0.39 
1 -,-1+0.24 
- L - 1J --0.21 


5 
4 


DWGl 
DWGO 


15 
15 


~ 
~ 




0.006 ± 0.002 
0.004 ± 0.001 


15.0 ± 2.3 
21.2 ± 2.9 


lmtl'l-l x io~ 3 

0.74 +1 n i 2 „ x 10~ 3 

— 0.53 


6.33+ 2 2 ;^ 


6 
6 



ROSAT 



ASCA 



20 
25 
18 
30 
30 
10 
13 



0.139 ± 0.042 
0.107 ± 0.035 
0.069 ± 0.040 
0.178 ± 0.037 
0.120 ± 0.050 
0.024 ± 0.013 
0.012 ± 0.006 



9.23 ± 0.57 
6.83 ± 0.73 
3.98 ± 0.23 
7.51 ± 0.72 
5.25 ± 0.76 
2.72 ± 1.37 
0.95 ± 0.23 



1.78 ± 0.21 
0.66 ± 0.02 
0.24 ± 0.02 
0.281 ± 0.042 
0.155 ± 0.037 
0.043 ± 0.009 
0.014 ± 0.005 



3.44 ± 1.81 
6.87 ± 3.68 
12.2 ± 7.3 
5.50 ± 0.96 
8.80 ± 1.63 
14.8 ± 2.5 
35.8 ± 10.2 



15.6 
5.83 
1.69 



+ 8.8 
-5.6 
' 2.25 



-0.75 

5.48 ± 0.62 
2.36 ± 0.31 
0.50 ± 0.13 
0.11 ± 0.03 



0.17 
0.75 
2.92 
0.65 
1.88 
6.23 
60.9 



+0.37 

0.12 
+ 2.16 

0.56 
+ 7.29 
-2.08 
+ 0.29 

0.23 
+ 0.97 

0.73 
+ 3.08 

2.35 
+ 57.9 

35.7 



References: (l)Carlberg et al. (1997); (2)Adami et al. (1998); (3)Mahdavi et al. (1999); (4)Borriello & Salucci (2001); (5)Sanders & Verheijen (1998); (6)van den Bosch & Swaters 
(2001); (7)Wu & Xue (2000); (8)Sato et al. (2000). 



fer significantly at low-mass end. The best-fit relations based 
on the X-ray sample alone are 

S c = l0 19Ai±2 - 65 (M vir /M Q r im±0 - 17 - : (1) 
c = 10 7 - 34±1 - 10 (M vir /M Q )-°- 45±0 - 07 . (2) 

In order to compare with theoretical predictions in hierar- 
chical clustering models, we compute the collapse redshift z co n 
for halo of mass M v j r identified at redshift z = in terms of 
the extended Press-Schechter formalism (Lacey & Cole 1993; 
Navarro et al. 1997): 

^crit (-^coll) ^crit (0) + 0.477 V2[<T 2 (/M vir ) - cr 2 (M vir )],(3) 

where / « 0.01, S cr [ t (z) = 8 cl - lt {0) / D(z, Om, ^a) is the den- 
sity threshold for spherical collapse at z with D being the linear 
grow factor and <5 cr it(0) « 1.69, and a(M v - l - c ) is the rms lin- 
ear density fluctuation at mass scale M v j r = (47r/3)p c riM-R 3 , 
which is related to the fluctuation power spectrum by 

<7 2 (M vir ) = — j dkk 2 P(k)W 2 (kR). (4) 

We use the standard top hat window function and parameter- 
ize the fluctuation power spectrum by P{k) — Ak n T 2 (k), in 
which the transfer function T(k) is taken from an adiabatic 
CDM model given by Bardeen et al. (1986) for the shape pa- 
rameter r = 0.25 and the Harrison-Zeldovich spectrum n = 1. 
Furthermore, we assume that the rms fluctuation amplitude 
within a sphere of R = 8 Mpc is <r% = 0.85. 

In terms of the prescription of Navarro et al. (1997), the 
characteristic density S c of a virialized dark matter halo M v ; r 
is proportional to the mean density of the universe at the corre- 
sponding collapse redshift z co n 

(M vir )] 3 , (5) 



where C is the proportionality constant and can be approxi- 
mately taken to be C = 3000 for our choice of cosmologi- 
cal parameters. In particular, for a power-law spectrum of pri- 
mordial density fluctuations, the above expression reduces to 
5 C oc M^™ +3 ^ 2 when combined with equation (3). There- 
fore, the best-fit <5 C -M v i r relation equation (1) from the X-ray 
data indicates n = —0.88 ±0.24, which is simply the combined 
result of Wu & Xue (2000) and Sato et al. (2000). The theoreti- 
cally predicted <5 C -M v ; r and c-M v ; r relations from equation (5) 
are illustrated in Fig.l. It is immediate that the data points of 
the galaxy sample (z < 0.3) are roughly consistent with the 
theoretical prediction. 

4. Discussion and conclusions 

Using the dynamical properties of dark halos reconstructed 
from galaxy distributions in massive systems (groups/clusters) 
and rotation curves in less massive systems (dwarf, low surface 
brightness and spiral galaxies) drawn from the literature, we 
have examined the 5 c -M v - lr relation predicted by hierarchical 
clustering models. It turns out that the observational data are 
in gross consistency with the theoretical prediction on a broad 
mass range from 1O 1O M0 to 1O 15 M , although the present 
analysis still suffers from the sparse data points especially on 
the low-mass scale of M < 1O 13 M . Indeed, this last point 
reflects the difficulty in the current determinations of the dark 
matter profiles of galaxies from their rotation curves. In fact, 
even if the high-quality rotation curves of HI disks can be 
traced out to large radii, the unknown mass-to-light ratio of the 
stellar components can result in ambiguity regarding the deter- 
minations of S c and r s in the NFW dark matter profile (e.g. van 
den Bosch et al. 2000). We have illustrated in Fig.l the influ- 
ence of M/L on the <5 C -M v i r and c-M v i r relations for DWG 
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sample. We cannot exclude the possibility that the best-fit val- 
ues of <5 C , r s and c in other two samples (UMA and DIG) may 
still contain rather large uncertainties due to the unknown M/L 
although M / L is allowed to vary in the fitting processes. Recall 
that good fits in DIG sample are obtained only for unreasonably 
large virial velocities and masses, while majority of the galax- 
ies in UMA sample simply failed in the fitting of NFW profile. 
So, more accurate and reliable data of galactic systems will be 
needed to tighten the constraints on the 8 c -M v - lr and c-M v j r re- 
lations on low-mass scales. Alternatively, the assumption that 
galaxies trace mass has been often used in groups and clusters 
in order to extract the dynamical properties of dark halos from 
the optically observed radial profiles of galaxies (e.g. CNOC, 
ENACS and PSG). A quantitative estimate of the uncertainties 
in the determinations of the NFW profile from this oversimpli- 
fication of galaxy distributions should be made in the future. 

The newly established <5 C -M v ; r relation from 'galaxies' dif- 
fers remarkably from previous findings based on X-ray mea- 
surements of the hot diffuse gas contained in galaxies, groups 
and clusters (Wu & Xue 2000; Sato et al. 2000). The latter re- 
sults in a shallower power-law spectrum (n w —1) than the 
expectations of typical CDM models on cluster scales. It is 
possible that the conflict is associated with nongravitational 
heating of the X-ray gas in the early phase of structure forma- 
tion. Such a speculation is supported by the good consistency 
of the 5 c -M v i T (or c-M v ; r ) relations between the two samples 
(galaxies and gas) in massive systems (M > 10 14 M Q ) and 
the remarkable difference on low-mass scales. Note that galax- 
ies are less affected by the possible existence of preheating. 
This situation is similar to the recent discovery of an entropy 
floor in the centers of groups and clusters (Ponman et al. 1999) 
due to the additional heating of the gas from supernovae and/or 
AGNs. Smaller systems (groups and galaxies) have lower en- 
tropies and hence are more affected than clusters. As a result, 
the present-day saturated configuration of hot/warm gas in less 
massive systems is shallower than the distributions of dark mat- 
ter and galaxies. Reconstruction of the dark matter profiles un- 
der the assumption of hydrostatic equilibrium without the de- 
tailed information about the gas distribution and temperature at 
large radii may contain large uncertainties. A quantitative anal- 
ysis will be needed to address whether the departure of the S c - 
M v ii relation detected in X-ray observations from those found 
by rotation curves and distributions of galaxies and predicted 
by typical CDM models arises from this speculation. 

Essentially, it is hoped that the <5 c -M V i r and c-M vlr rela- 
tions can be extended to less massive systems below 10 10 A/q. 
This last point is potentially important for distinction between 
CDM and warm dark matter models (Eke et al. 2000). Alter- 
natively, the employment of the 8 C -M V - U relation established 
in the present work can allow one to test the different prescrip- 
tions of halo formation time (e.g. Navarro et al. 1997; Salvador- 
Sole et al. 1998; Eke et al. 2000; Bullock et al. 2001), and we 
will report the result elsewhere. 
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Fig. 1. The observationally determined <5 C -M v i r (upper panel) 
and c-M v ; r (lower panel) relations for two samples: (^galax- 
ies - CNOC (open triangle), ENACS (open circle), PSG (open 
cross), DWGO, DWG1, DIG and UMA (open squares from left 
to right); (2)X-ray gas - ROSAT (fitted square) and ASCA 
(filled circle). The asterisks from left to right denote the re- 
sults for five individual objects: DDO 154 from the measure- 
ment of HI rotation curve (Carignan & Purton 1998; Burkert & 
Silk 1999), the Hydra A cluster observed recently with Chan- 
dra (David et al. 2000), the poor cluster AWM 7 from the NFW 
fit to the radial profiles of galaxies and their velocity dispersion 
(Koranyi & Geller 2000), the nearest cluster Virgo from a com- 
bined analysis of optical and X-ray observations (McLaugh- 
lin 1998), and the cluster Abell 576 from the measurement of 
galaxy distribution out to the infall region (Rines et al. 2000). 
Dashed line is the best-fit relation to the X-ray data. Solid line 
is the theoretical expectation for a flat cosmological model of 
n M = 0.3 and Q, A = 0.7 with T = 0.25, cr 8 = 0.85 and n = 1. 



